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Abstract 

The new symmetry of the loop dynamics of QCD is found. This is a local SUSY 6s = 
9f3{s),50 = j3{s) of the superloop field X^{s,0) = x^{s) + 97p^{s). The remarkable thing is, 
there is no einbein-gravitino on this theory, which makes it a ID topological supergravity, or 
locally SUSY quantum mechanics. Using this symmetry, we derive the large Nc loop equation 
in momentum superloop space. Introducing as before the position operator we argue 
that the superloop equation is equivalent to invariance of correlation functions of products 
of these operators with respect to certain quadrilinear transformation. As a consequence of 
this nonlinear symmetry, the coefficients of the Voiculesku expansion of the position operator 
satisfy recurrent equations. The generators of this symmetry are involved in the glueball 
spectrum, as it follows from the loop-loop correlation function equation. The correlation 
functions of external flavor currents in the background of constant flavor gauge fleld with 
flnite density of topological charge in the chiral limit are also considered. We represent them 
as certain flnite dimensional integrals in superspace, involving the Greens function, which is 
expressed in terms of X^. The expansion in powers of external momenta A;i . . . A;^ is calculable 
in terms of the above universal numbers in the Voiculesku expansion. 
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1 Introduction 



Recently |[T[, some progress was achieved in a long-standing problem of string representation 
of large N QCD. We introduced the position operator of the endpoint of the QCD string 
and applied the momentum loop equation (MLE) to get some relations between the 

planar connected moments of vacuum expectation values of this operator. 

The non-commutative probability theory j^, |5|, |^] was then used to build the Fock space 
and operator expansion of in powers of the creation operator aj^ which obeys the Cuntz 
algebra a^aj, = 6^,y. The Fock space consists of words aj^^^aj^^ . . . aj^^ |0), where the vacuum |0) 
is annihilated by a^. 

There are several problems with this approach. One problem is that some coefficients of 
this operator expansion 

= + ^«]^ + Balalai + C (4aJ,aj, + aj,aj,4) + . . . . (1) 

are left undetermined. Apparently, some information is still missing. Another problem is that 
there seems to be no direct relation between these coefficients and the QCD observables. 

As for the first problem, it was known since the first papers on the loop equation |3| that 
one should restrict oneself to the class of so-called Stokes type functionals to reconstruct the 
Wilson loop from the loop equation. In those papers it was shown that iterations of the loop 
equation within this class of functionals reproduce the Faddeev-Popov perturbation theory 
including ghost loops. 

The Stokes type functionals by definition satisfy the Mandelstam relation 

^ ^U^)ir^> (2) 



he area derivative ^ ^ ( ^ is ^ 

identity 



where the area derivative ^ is a skew symmetric tensor operator satisfying the Bianchi 



d 

•5" (•5) ^ TT + cychc = 0. (3) 

The derivative is defined as 

/•s+O ^ 

W= dt—— (4) 

Js-o dx^[t) 

These relations were not taken into account in our previous paper Also, the loop 
equation used in that paper was not the most general one. It corresponded to the original 
vector loop equation |Q] multiplied by the tangent vector x'^{s) and integrated around the loop. 
Normal component terms were eliminated by this projection. For example, in two dimensions, 
the loop averages were computed Q using the normal component of the loop equation. The 
tangent components in this case were identically satisfied for the area dependent Ansatz. 

The second problem is more serious. The operator X as defined is apparently not an 
observable, so it may be infinite or zero in the local limit of the theory. 
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In this paper we fix both of these problems at the same time. We choose the superloop 
formalism (see, e.g. ^]) which allows one to impose Mandelstam and Bianchi relations 
as local supersymmetry relations. We derive the large superloop equation in momentum 
space. 

The SuperFourier transformation we use here, is free of the zero point fluctuation problems, 
which were a serious obstacle to the old momentum loop dynamics. These fluctuations cancel 
by themselves, without any Wiener measure, by virtue of supersymmetry. 

Thus, there are no more hidden infinite renormalization factors in our momentum super- 
loops. The corresponding super symmetric generalization of the position operator of the string 
obeys a different set of recurrent equations, which are more restrictive than before. 

These relations are based on the important observation that there is a local SUSY in 
the loop dynamics. The local SUSY is valid in the large Nc loop equations as well as in 
the relations for these observables in the chiral limit. We are able to explicitly find all the 
superspace-dependent quantities (they reduce to the SUSY B functions). 

The physical interpretation of these relations is found. They represent the symmetry 
relations for an operator X correlation functions. One can perform certain nonlinear trans- 
formation of X operator with two vector parameters. This symmetry is completely equivalent 
to the momentum loop equation. 

The non-commuting vector generators of this symmetry act as "QCD string" endpoint 
translations. Their eigenfunctions must have fixed spatial momentum, which condition yields 
the glueball spectrum equation. We derive this condition from the glueball superloop equation. 

The superloop in momentum space has direct physical meaning in the chiral limit of QCD. 
It is related to the correlation functions of external flavor currents. The local supermomentum 
is expressed in terms of the set of external momenta entering the quark loop. 

One set of correlation functions involves the chiral anomaly (the periodic boundary con- 
ditions for the Fermi component of the superpath). Those are expected to be completely 
dominated by the perturbative contribution. The argument is based on the index theorem. 
However, in presence of non-smooth fluctuations of the gluon field this perturbative argument 
may fail. We relate the coefficients in the low energy expansion of these correlation functions 
to the coefficients of the Voiculesku expansion. 

Another set of calculable observables is given by the ordinary correlation functions of ex- 
ternal flavor currents taken in the background of external field with finite density of topological 
charge (the anti-periodic boundary conditions for the Fermi component of the superpath). The 
coefficients of the low energy expansion are also related to the coefficients of the Voiculesku 
expansion. We derive an explicit set of operator equations in presence of external field. 

The resulting theory can be reformulated as the theory of the quark superfield propagating 
around the superloop, rather than target space. The corresponding propagator can be found 
exactly in terms of the quark position operator X^. This resembles the string theory, but we 
do not see any precise correspondence. 
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In order to handle the chiral symmetry breaking in absence of the external density of 
topological charge we included the quark mass term. This resulted in extra term in the quark 
effective action on a superloop. This term appears to be a quartic interaction, rather that a 
mass term. This superfield theory is essentially one dimensional, so it will hopefully present 
no new problems. 

2 Vacuum Energy and Superloops 

Let us present our superloop formulation of QCD. It starts with the formula (derived in 
Appendix A) for the quark loop energy in the presence of color and flavor fields 

£ = const - ^tr log (mo + (iT^iV^,)^) = (5) 
^ dTT"^-^ (1 + alnT) e-'^<ti exp (-T {i-f^V^f) . 

In this and subsequent formulas the limit a +0 is understood. This a prescription is 
equivalent to the C, rcgularization and it goes along with dimensional rcgularization. 

The covariant derivative operator involves all the gauge potentials V ^ = + + B^, 
where Afj_{x) is the quantum gluon field, and B^(x) is the external flavor fleld. 

The SUSY follows from the effective Hamiltonian 

itl^Vi,f = {QL,QR}, (6) 
Ql,r = («7mVm) 

Therefore a Dirac operator generates the (global) SUSY transformations. This is the well- 
known SUSY Quantum Mechanics. 

The path integral representation can be formulated in terms of the super trajectory 
a;ju(s), V'//(s). The Fermi coordinates V'/i(s) represent the the matrices 7^\/T. 

It is most convenient to use the superspace S = {s,6) with the superfleld 

X^{S)=x^{s) + 9ij^{s). (7) 
with boundary conditions of (anti)periodicity 

Xf,{l)=x^{0),^P^{l) = -^P^{0). (8) 

The negative sign for V'/i is needed to make the traces of 7 matrices cyclic symmetric in spite 
of anti-commutation of -i/;^!. Moving the left to the right all the way through remaining 
(odd) number of ip would produce the minus sign, unless we compensate it by the extra sign 
flip at s = 1. 

The action is given by an integral over superloop 



\s J d9C{S) (9) 
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with the super Lagrangian (in our normahzation of proper time and ip fields) 

C{S) = £'=^"(5) © DX,,{S)A^{X{S)) © DX^{S)B^{X{S)), (10) 
C'^-iS) = -j^DX,{S)D'X,{S) 

Here 

d d 

is a covariant derivative in superspace. Its square is the time derivative 

= (12) 

OS 

The direct sum of the color and flavor matrices in understood in the sense of taking the 
ordered exponential of the action in the (super)path integral 

J PXrexp(J) = J PXexp (^J^ ds J dO C^'^'iS)^ (13) 

Texp (^j\s J dODX^iS) {A^{X{S)) ® B^{X{S)) 

which is the non-Abelian generalization of the Feynman's path integral. One may prove that 
the results are equivalent to the more traditional approach with the ordinary exponential of a 
second quantized action. 

Conventional formalism shifts the momentum by =^ zp^ © Afj_ © after which the 
gauge fields go to kinetic terms in the Lagrangian, which is inconvenient. So we undid this 
shift, coming back to original ideas of Feynman, who started with the (Abelian) Wilson loop 
as an amplitude for electron propagation. In our previous paper Q we discussed these issues 
in detail. 

3 Superloop Kinematics 

In our case it is implied that the gauge fields have the superfield as their argument, which 
means that first derivative terms are also included. The (global) SUSY generalization of the 
the ordered product is also known |]l^]. One introduces a SUSY generalization of the distance 

L{S, S') = s' -s + 96', (14) 

and the corresponding theta function 

e{s, s') = ems, s')) = e(s' - s) + 9e'6{s' - s). (15) 

The covariant derivative D = -§0+9-^ of this Q function reduces to the SUSY 6 function 

DQ{S, S') = {9' - 9)6{s' -s) = 6{S' - S). (16) 
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This S function is a Fermi rather than Bose element, and it is an odd function 

5(S) = ~6{-S), (17) 
S{0) = 0. 

The last property is very convenient for computations of the SUSY graphs. 
The generalized function satisfies the usual property 

Q{A,B) + Q{B,A) = 1, (18) 

from which it follows that @{A,A) = i. 

We shall also need the generalization of the finite range integral 

f-B roc r 

/ dS= ds dee{A,S)e{S,B), (19) 
as well as the multiple ordered integral 

/ = / dsi d0i... dSn dOnQiA, 5i)e(5i, S2)... @iSn-l, 5„)e(S„, B). (20) 

J A J —oo J J —oo J 

The following superfield differential forms will play an important role 

dX^{S) =dsdeDX^{S), (21) 

dJ^X{^y{{S})=dX^^{S^)...dX,^{Sn). 

By construction these forms are globally SUSY, later we prove that they are locally SUSY as 
well. We are using the curly brackets {S} for lists of indexes or arguments and the tensor 
power notation 

n 

^w({^})^n^M,(^0- (22) 
1=1 

The superfield differential 1-form satisfies the usual integration identities 

f'' dFiS)=FiB)-FiA), (23) 
J A 

d{F{S)G{S)) = dF{S)G{S) + F{S)dG{S), 
which allow integration by parts. The finite range integrals are additive 

f-B fC f-B 

/ dF{S)G{S) = / dF{S)G{S) + / dF{S)G{S). (24) 
Ja Ja Jc 

The integrals of multiple forms satisfy the composition law 

r dp,{s) [' tf"P{^}({c/}) r <rP{.}m) = r d^+^+'P{,}{{w}), (25) 

with the joined list of indexes 

{p} = {M,A,M}. (26) 
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as well as the multiplication laws, the simplest of which is 

^ciPxiU) r d^P{,}{{V}) = X: r cr+'P{.(i)}{{S}l (27) 
J A J A 1—qJA 

where 

{^(0} = {/"!>•• (28) 

In other words the index A is inserted at the position Z + 1 in the list of /x indexes. 

Both of these laws follow from the local SUSY which we prove in the next Section and the 

supersymmetrization theorem which we prove in Appendix C. 

Using these forms, we can define the SUSY path exponential with initial point A = (a, 6a) 
and final point B = (b, Of,) 

U^[A„X,]^TeKpi dX^{S)A^iX{S)))^J2 ({^(^)}) ' (^9) 

One may verify the multiplicativity of the path exponent 

[A^, X^] = U^[A^, X^pi [A^, X^] , (30) 
its differential with respect to the endpoints 

dUi[A^,X^] = 5{S - A) +Ui[A^,X^]dX^{S)A^{X{S)), (31) 
dUi[A^,X^] = -6{S -B)- dX^{S)A^{X{S))Ul[A^,X^], 

as well as the usual formula for its variation with respect to the gauge potential 

5a,K[A^,X^] = r dX^{S)Ui[A^,,X^,\5A^{X{S))Ui[A^,,X^,\. (32) 

J A 

In particular, for the gauge variation 5gaugeA^{X) = V^q;(X) (note that the gauge parameter 
a{X) having the superfield -'^(>S') as its argument) 

5gaugeUf[A^,,X^\= ^ dX^{S)U^A[A^, X^]V ^a{X{S))U§ [A^, X^]. (33) 

J A 

Using above endpoint derivatives we reduce the integral to the total derivative, up to the 5 
terms, which yields 

5„Wf = Ui[A^,XMX{B)) - a{X{A))U^[A^,X^]. (34) 

We got something new here. The trace of this variation will vanish provided the superloop 
is closed, i.e. X{B) = X{A). The Bose and Fermi parts by themselves could be opened. The 
gap in Bose part x{b) — x{a) = —9h^{b) + Oai^ia) represent the nilpotcnt even element of the 
Grassmann algebra. Its higher tensor products {x{b) — x(a))^^ . . . (x(6) — x(a))^^ vanish at 
n > 2. This is less restrictive than the ordinary gauge invariance which would demand n > 0. 
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We are in a position to define the superloop 

^ = ^ (trZ^|[^^,X^] . (35) 
In usual notation the SWLoop reads 

W [X]l = ^trf exp ds (^A^{x)x'^ + ^^^V'.F^,(x))) ^ W [x, V^] • (36) 
This formula was first obtained in [12, 13 1 for the same purpose of supersymmetric description 



of spin I particle. In case of the Grassmann variables at the end points A,B, there will be 
extra terms. Say, in Abelian case there will be extra factors exp (ib0^/'^^^(x)) at the endpoints. 

The relation between this formula and the one with superfield is quite amazing. The 
ipipdA part of the il^tpF terms comes from the expansion of the potential in the action, and 
the commutator terms ipipAA come from the 6i9j6{si — sj) term in the SUSY definition of 
the the ordered product ||ll[]. The SUSY and the non- Abelian gauge invariance work in close 
collaboration, which is quite common, as we shall see. 

The vacuum energy can be written in terms of superfields as follows 



£ = -^1 (iTr"-^(l + alnT)e-^""o (37) 







JvXexp(^-^J^ dX^{S)D^X^{S)yvfexp(^J^ dX^{S)B^{X{S))^ W [X]^ . 

The ordinary perturbation theory corresponds to expansion in both flavor and color gauge 
fields and taking the Gaussian path integrals term by term. The role of the Grassmann part 
G'lpnis) of the path X^(S') is to get correct numerators of the Feynman diagrams, coming from 
the Dirac traces in usual language. 

This follows from the explicit form of the kinetic energy 

^ J^' dX^{S)D'X^{S) = ^ ds {x[,{s)x[,{s) + i''^{s)Ms)) ■ (38) 

One could rescale s ^ T s^ip^ ^ ^/TtI}^, which only affects this term, as the rest of the terms 
are parametric invariant. Then the kinetic energy will take a familiar form, which corresponds 
to the usual free particle. The il) part after quantization reproduces the usual Clifford algebra 
of Dirac matrices. The extra factors of mass in the local measure, coming from rescaling of 
V'/i, are just the correct factors needed to normalize the Wiener measure for x^{s). 

We mentioned this only to make our path integral look more familiar. For our purposes 
this transformation is not needed, in fact, it would only make the formulas look more complex 
than they actually are. 



For example, the Euler-Heisenberg effective Lagrangian |14] in homegeneous external field 
By{x) = I B^yXfjt is computed in just few lines in Appendix B. 
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4 Correlation Functions of Vector Currents 



Let us expand the vacuum energy in powers of the external field. The expansion coefficients 
are related to the correlation functions of the vector flavour currents. 
We use the Fourier integrals 

f df^k 

B>.{X) = J j^B^{k)e^p{zk^X^), (39) 

which also define the SUSY extension to X = X{S). In the n-th order we shall have the 
product of n exponentials, which can be written in terms of the supermomentum 

exp (^i ^ kiX^iSi)^ = exp £ dK^{S)X^{S)^ = exp (^-i dX^{S)K^{S)j , (40) 

with the supermomentum 

K^{S) = J2ki,@{Si,S), (41) 

i 

dK^{S) = Y^ki,dSd{S -Si). 

i 

We need the standard SUSY functional derivative for a superfield Pn^s) = Pij,{s) + 6ip^{s) 

' ' +0A^, (42) 



6P^{S) Sip^{s) 6p^{s)' 
which satisfies the identity 

^J^^P,(S')=6,,S{S-S'). (43) 

In the same way, as we introduced the superfield 1-form dXn{S) we could introduce an invari- 
ant superfield gradient 1-form 

^ dSTT^- (44) 



The motivation for this notation is the following identity 

d 



dP,{S) 



f dX^{U)P^{U) = dX^{S), (45) 

J A 



which means that this gradient form acts on the integral of differential form just as it does in 
Bosonic case. 

The vacuum energy takes the form 

£ = tfj^,---l^,^' (j:k}iB;iik,)...BZiK)(j;il...JZ)ih...K), 

(46) 
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with the following planar connected correlation functions of flavour currents 



j;j...J2:)(/ci...M= (47) 

dM [K]l 



,iS)J Jo dK^,{Si)...dK^„{Sn 
where we introduced the following operator function 

Z{L) = lim n dTT"-^ (1 + a In T) exp ( -Tm?. - - | . (48) 

a-»0 Jo V / 

This is a certain Bessel function of its operator argument. 
The loop functional 

M [P]^ = J VXS'^iXiA)) exp dP^{S)X^{S)^ W [X]^ ; (49) 

is a Fourier transform in superspace. 



5 Chiral Anomaly 

There is a simpler observable, which is not influenced by quark kinetic energy. This is the 
chiral anomaly 

I[B] = { tr 75 exp (-T (z T^V^)') ) • (50) 

as a functional of the external flavour field B. The famous index theorem states that it is 
T-independent: due to the 75 symmetry of the Dirac operator, all the finite eigenvalues drop 
from the trace, so that only the T- independent contribution of the zero eigenvalies remain. 
In absense of quantum fields, this allows to compute it, taking the limit T ^ and using the 
WKB expansion. This leads to the the Grassmann integral over the zero mode 

Ici[B] cxvj dVtrTexp Q5^.V|!V'°) oc Ve^,...^,ti B^,^, . . . B^^_^^^. (51) 

In case of abelian field this is the pfaffian VdetiB. 

The quantum chiral anomaly, I[B] according to the same index theorem, is simply the 
vacuum average of the classical chiral anomaly of the direct sum of the two gauge fields 
( Ici[B © A] ) computed in absence of the B field. 

I[B] = {Ici[B®A]) = NJa [B] + Nf { la [A] ) . (52) 

The last term vanishes by parity, so that the chiral anomaly stays the same as in the free 
quark theory. 

Let us now relate this vacuum average to our momentum loops. We have to repeat the 
steps of the previous section, with two changes. First, we choose the different boundary 
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conditions. Now both x{s) and iJj{s) are periodic, to compensate for the 75. After that we can 
drop the kinetic energy term, according to the index theorem. 
We find the following formula 

NJcl[B] = ^ / PXfexp ( [' dXi,{S)Bi,{X{S))) W [X] , (53) 

J periodic \J0 / 

which provides an important normalization condition for the SMloop considered below. 

In particular, consider the case of constant non-abelian field Bijl{X) = const. This corre- 
sponds to constant field strength B^,^ = [B^ , Bi,\ . The chiral anomaly reads 

Bfx^ B^^ . . . B^^ . (54) 

Comparing this with above path integral we get the sum rule 

/ VX [' d%^y {{S}) W [X] oc (55) 

Jperiodic JO 

This fixes the coefficients in front of Kf^^ . . . K^^ in Fourier transform of W [X\ with periodic 
boundary conditions. 

6 Local SUSY 

Let us study the symmetries of the SWloop. It is invariant under parametric transformations 

5x^{s) = a{s)x'^{s), d%l)^{s) = a(s)V';[j(s) + ]^Oi {s)^} ^,{s) , (56) 
which is straightforward to verify. In addition it is invariant under local SUSY transformations 

5x^{s) = (5{s)Ms), 6^^{s) = P{s)x'^{s). (57) 

with some Grassmann variable P{s). Together these transformations define the reparametriza- 
tion of the superfield, 

5X1, = Xf^is + 6s,e + 59) - X^{s, e) = (^6s^ + Se^^ X^{s, d), (58) 

with 

Ss = a{s) + ep{s), 50 = ]^a.\s)e + P{s). (59) 

Prom now on, we shall always mean the local rather than the global SUSY. In Abelian 
theory this is just the formal symmetry of the gauge field part of the action. One may verify 
that the variation of the A— term in a Hamiltonian after integration by parts cancels the 
variation of the F— term. In the last variation the extra terms with gradients of F^i^ add up 
to zero by virtue of the Bianchi identity 

^i.i^AxdxFi,^ = (60) 
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In the non-Abelian theory the same can be proven for the the ordered exponent. In the 
superfield language the SUSY of the action follows from the opposite transformation laws of 
the integration measure and covariant derivative 



6dS 

and 



d(Ss,5e) \ ,^ /d5s d69\ ,^ fl , A ,^ 

so that the differential d = dSD stays invariant. 

It remains to check the function (|l5|). The parametric invariance is obvious, as for the 
local SUSY, it can be verified as follows 

5G(5, S') = 6L{S, S')6{s - s' + 66') = (63) 
{6 + 6'){(3{s') - Pis))6{s -s' + 66') = 
{6 + 6'){(3{s')- l3{s))5{s-s') = Q. 

It is interesting that Andreev and Tseytlin [|^], who introduced this SUSY G function, did not 
notice its local SUSY, as they were working in a superstring theory, where this local symmetry 
was broken down to a global one. 

From the point of view of the ID supergravity this local SUSY corresponds to the fact 
that the Wilson loop does not depend on the metric in ID superspace, described by the 
einbein-gravitino superfield. Thus it measures topology in this space, which is rather simple. 

It is quite instructive to verify the local SUSY directly, using the components of the 
superfield. We refer the reader to the review paper where the loop kinematics and dynamics 
was discussed in great detail. Using the methods of that paper one can prove the SUSY of 
the non-Abelian SWLoop as follows. 

From the point of view of loop dynamics the SWLoop is not an independent functional. 
It relates to the Wilson loop by the linear operator 

W[x,V']o = rexp ^'dsiv^V'.^^ W[x]l. (64) 
Using this representation, the functional derivatives can be readily computed 

(65) 

^ ^ (4(»)^». + !.„(»)*(»).,.(.)) ^(.). 

It is implied that the left functional derivative is used for the Grassmann variable ip. 

Multiplying the last relation by ip^is) and using the non-Abelian Bianchi identity we 



dXfj,{s) " ^ da^i, dtp^ 



(66) 
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which is a differential form of the SUSY relation. 

This local supersymmetry is a genuine symmetry of the loop dynamics in the pseudoscalar 
sector. It does not require any external einbein-gravitino fields, which were needed to super- 
symmetrize the massive spin ^ particle. In other words, we are dealing with topological SUSY 
Quantum Mechanics. 



7 SuperLoop Equation 

The SWLoop equation is obtained as a SUSY extension^ of the usual loop equation |||, |3|. It 
looks exactly the same as the ordinary loop equation, with a superfield Xfj,{S) in place of the 
ordinary coordinate x^{s), 

° = r [jJ^M^BiS) ^ DX,{0)DX,{S)5''{X{S) - X(0))W [xf, W [X]^^ . 

(67) 

In Appendix C we derive the supersymmetrization theorem, which states that any locally 
SUSY functional of the curve T : S = {s{t),9{t)) which vanishes identically for F being the 
real axis: 9{t) = must vanish identically in the whole superspace. The bottom line is, the 
local SUSY transformation has one Bose and one Fermi function of one variable, so it allows 
us to move and reparametrize any curve from the real axis anywhere in superspace. 

This resembles the analytic continuation from the real axis to the complex plane. As long 
as the Cauchy-Riemann equations are satisfied this analytic continuation looks like a trivial 
replacement of the real number by a complex number in any equation. The hidden cost is just 
the Cauchy-Riemann equations which should be checked in every case. 

One may think that the SUSY loop equation follows from the ordinary loop equation by 
means of this naive SUSY extension, but this is not true. The supersymmetrization theorem 
does not literally apply here. The two arbitrary functions of the SUSY transformation are not 
sufficient to eliminate the Grassmann components of d- vector superfield X^{S). One has to 
rederive the equation from scratch, using the superfields instead of coordinates. 

Let us examine the linear term in components 

' aT^TH^?^ = -^/^(O) a7^^ [^lo = -V'.(0)V^F^.(x(0)) W [X]l . (68) 

with the colour matrix V^F^jy(a;(0)) being inserted in the beginning of the the ordered product. 
According to the quantum YM equations, this can be replaced by 

with the same agreement about insertion inside the the ordered product. The variation of the 
path the ordered exponential goes as follows (each step can be proven using the above Taylor 
'^I am indebted to Sasha Polyakov for the advice to switch to superfields here. 
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expansion with SUSY @ functions) 




a 



1 




(70) 



{S)S\X{S)-x{0)) 



Texp 




The function here is the usual delta function in d dimensional Euclidean space, with argu- 
ment replaced by the superfield. It should be understood in the sense of a Fourier integral 



Recalling that ^',^(0) = DXi,{0) and Xi,{0) = X^{0) we get the above superloop equation. 

Let us see how this works in components. The integral over the 9 variable in (iX^(S') 
picks up the first derivative in the superfield expansion. This is equivalent to applying the op- 
erator V^(s)5^(s) to the d'^{X{S)-X{0)) function as well as to the SWLoops W [X]^ W [X]^. 
In the latter case these are the covariant derivatives d^{s). 

Like in a SWLoop itself, the various parts of the covariant derivative come from various 
places, supported by local SUSY. The ^ terms in covariant derivatives come from expansion 
of the X superfield, and the terms come from the expansions of the G functions in the the 
ordered product (the last @ function @{Sn,S) in the W [X]^^ and the first one @{S,Si) in 
the W [X] are differentiated). We get simply 



as an operator applied to the rest of the factors. 

The same expression could be obtained by an honest variation of the SWloop. In superfields 
everything looks so simple, but in fact all the subtleties of the loop dynamics are automatically 
taken care of. 

Let us discuss the boundary conditions. The initial loop W [X]q was periodic in x{s) but 
the boundary values of 'ip{s) were left free. The loops W [^]o )^ i-^]^ nonlinear part 

of equation belong to the same class. 

The right boundary value ip{s) of the of the first part of the path coinsides with the left 
boundary value for the second part, so that we have a convolution in the loop equation. 




/d k f 
J d'^q5'^{k - q) exp (i A:^X^(5)) exp {-i g^x^(O)) 

— ^ exp (z k,{X,{S) - x,m) = 6\X{S) - x(0)). 



(71) 



(72) 
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The physical loops, which enter the vacuum energy and the vacuum topological charge, 
are closed in ip with antiperiodic or periodic boundary conditions. This can be achieved by 
setting = ±^(0) and integrating over ^{0). 



8 Momentum Superloops 

Let us now go to supermomentum space. There is one important point to clarify before we 
do it. The natural definition would imply the periodic momentum P^{S) , P^{1) = Pfi{0). 
However, the loop equation does not close in the space of periodic momentum loops. 

This issue was discussed before ^, ||. In general, one could introduce some gaps AP{Si) = 
P{Si + 0) — P{Si — 0), in which case the momentum loop equation relates the n— gap functional 
to bilinear superposition of / + 1 and n — I + 1 functionals, with / = 0, . . . n. 

These gaps correspond to slow decrease ^ of the Fourier harmonics of p^{s), which is an 
unpleasant complication of momemtum loop dynamics. As we see it now, this complication 
can be avoided. Take non-periodic functions, -P/^(l) 7^ -f/x(0) which are otherwise smooth. The 
derivative DP^{S) could be finite everywhere including the endpoints. However, the endpoint 
derivatives are different I?P^(0+) / DP^{1—). 

From the point of view of the general framework, with gaps, this corresponds to only one 
gap, between 1 and 0. The inside momentum loops will also have one gap each, between their 
initial and finite points respectively S, and 1,5". In other words, the "loop" has a topology of 
an interval, rather than a circle. Cutting the interval at some inside point we get two intervals, 
so that the topology is trivially preserved. 

The Fourier transformation is defined as follows 

M [P,V',V'']J = / VXexp (^i j\^{S)dP^{S)^ W [X] , (73) 

where ip,ip' are boundary values of the Fermi part of the path X^(5'). 
After the Fourier transformation the loop equation reads 

'li /^^^.c^_ r,^ I ,,,,.dM[P,^,^{s)f,M[PMs),^fs 



DP^{o)M jl<^P^^s) = -^-al f^dS J dU{s 



6P^{0)6P^iS) 



(74) 

This equation is not as simple as it looks. Several things needs too clarified. 

The coordinate 5 function disappeared in the same way as in the ordinary momentum loop 
equation §. The product 6'^{X{0))6'^{X{S) - X(0)) can be replaced by 6'^(X{0))6'^{X{S)) 
which is just what is needed to get the correct integrand in the product of the two Fourier 
transforms. The 6'^{X{0)) function eliminates translations in W [-^]o the other one, 
6'^{X{S)) does the same for M[X]l. 

Note that in the last delta function the 9 term in superfield is present, which adds the 
terms with gradients of the delta function. This is not a problem, since the same SUSY 6 
function was included in the definition of the Fourier transform. We had to do it to preserve 
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the local SUSY. By virtue of SUSY the choice of the point to fix on a loop is arbitrary, so 
that we could replace d'''{X{S)) by 6^{X{l))m the Fourier integral for M [P]s- 

Some subtleties of the factorization of the Fourier integrals are discussed in Appendix D. 

Note, that there was no delta function for the Fermi part of the loop, therefore the inte- 
gration over ^p{s) was not eliminated. There is no explicit dependence of these values in our 
equation, just the integration. Therefore, we can look for the diagonal solution 

/I \ /I \ 



M [P, ^, i^TA = ^ [2^^ + ^')) ^ [P+]a + [-^(^ - V'O j M [P-]^ . (75) 
Here M correspond to the two types of boundary conditions 

M [P±]^ = J d^i^M [P, V, ±i^]A ■ (76) 

We used the fact that (5'^(0) = for the Fermi variables. 

Thus M. [P+\ describes the chiral anomaly, and M. [P—\ - the vacuum energy (in the latter 
case one should also include the quark kinetic energy - see below). 

Then, the delta functions integrate out, and we are left with the single equation 

np.m [Pli £ ,PAS) = 2N^,i [ -isSS^, (77) 



where we introduced the chiral functional 



C[P] = \{M[P-]+M[P+]) (78) 

which corresponds to projection operator | (1 + 75) inserted into the Dirac matrix trace. This 
projection operator commutes with the Lorentz generators 

(y^v = ^ [7m > lA > (79) 

which enter the Pauli term a^nF^y in the square of the Dirac operator. Hence, the chiral 
functional corresponds to the chiral generators ^ (1 + 75) u^,^. The usual even/odd functional 
correspond to the even/odd parity part of the chiral functional. 

In spite of chiral loop equation being parity even, we have to include the terms in 
the coefficients. These terms are fixed by the chiral anomaly, as discussed above. We shall 
come back to this issue below, after we introduce external field. 

The SMLoop equation can be simplified if we introduce the SUSY area derivative 

or, in terms of differential forms 

9C[PfA _,^(^. dC[P]^^ 
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Later we compute the area derivatives for the terms of the SUSY Taylor expansion of the 
SMLoop. As discussed before, we assume the smooth non-periodic function P{U), so that 
P^{S+) = P^{S-) but P^(l) / P^(0). 

In the SMLoop equation both hnear and bihnear terms involve the (iP^(O) = ) + 

^dPf^{0+) form. Both of these forms being arbitrary, we could "cancel" them, i.e. leave only 
the tensor coefficients multiplying each form. 

This yields two equations 

where ^ = or ^ = 1. 



9 The External Field 

The problem with the momentum loop dynamics is that it is non-perturbative. The zeroth 
approximation of perturbation theory corresponds to singular momentum loop which singu- 
larity we expect to disappear in the full solution. The only way we can study this full solution 
so far is to expand in powers of momenta. The expansion coefficients are expected to be 
proportional to the powers of confinement radius which is infinite in perturbation theory. 

The loop equation being nonlinear, we can never know whether the expansion we build, 
really matches the perturbative QCD as it should at confinement scale of momenta. What if 
this is a wrong solution, some strong coupling artifact, like the ones in lattice gauge theory? 

There is a following way around this obstacle. Add the external Abelian constant field 
Bfj^iy, then at large B^i, the asymptotic freedom will hold, and we will know the solution. At 
any finite field the SMLoop will be still expandable in powers of momenta. 

The B- dependence of expansion coefficients can be found from the following identities 

f^ft<iS.e(S.5.),^#M^. (83) 



These identities follow directly from the original functional integral, with the Abelian external 
field factor inserted 

fexp (^i dX^{S)B^{X{S))j = exp (^^ B^, Q{Si, S2)dX^{Si)dX,{S2)^ . (84) 

In the Fourier integral one could replace dX^{S) by i sp^{s) • 

It is convenient to introduce the area derivatives, then one finds in (^) after some algebra 



dC[P]^ 



A _ M ^ , ^ \r [Die 



dB„ 




C[Pn+ (85) 



' dPx{Si)dPs{S2)e{Si,S2) ^ 
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As for the loop equations, they are modified in an obvious way, by adding the derivatives 



of this exponential to the operators 



5 



The resulting equations read 

= dP^iS) ( ( 5pA + * Bp^-^ ) Up^ + ^ Bps^ 1 C [P]' + (86) 



2 dC[PY^c[P]l 



where ^ = or ^ = 1. 

It is convenient to pull out of the SMLoop the exponential, corresponding to the free quark 
theory 

C [P]^ exp (^1 P,{S)dP,{S)^ C [P]^ . (87) 

This factor being multiplicative around the loop, it will cancel on both sides of the loop 
equation. The only result will be the shift 

This shift eliminates the DP term in the linear part of the momentum loop equation, and 
we get 

B,^B,s C dP^{S)-^-^C [P]l = (89) 

N.d /; ^P.(S) (. + ^) (. B.-' + s^) C C [P| J . 

Now, the perturbative solution would correspond to iterations of this equation in gQ, start- 
ing with the free quark loop. The free quark computations are presented in Appendix C (see 
also below). The result for initial definition of chiral supermomentum loop functional reads 

C [P]^ ^ VdetiB 9 (VdetzS) exp (^^ B^^ Pp{S)dP,{S)^ + O [g^) . (90) 

It is an interesting challenge to verify how the higher terms of perturbation theory in a 
background field are reproduced. 

10 The Position Operator and Nonlinear Symmetry 

Let us proceed with the momentum superloops. The next step in Q was the introduction of 
the position operator as a connection in momentum space. 
The corresponding generalization of this Ansatz would be 



f exp L j\p^{S)Xp 



) . (91) 
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The operator X„ is an ordinary Bose operator. There are no restrictions on the commu- 



tators 



X^ , Xy 



so that these operators form a free algebra. 
This Ansatz involves the superloop in momentum space with constant non-Abelian con- 
nection iX. In components: 



Texp [ z I dP^{S)Xf}\ =fexp(f ds (iX^p'^ - ]- 



Xfj, , Xi, 



(92) 



The above definition of the position operator X does not assume periodicity of the superloop. 
One could translate X by a constant vector operator C^, which commutes with X^ as well as 
itself: 



X,. 



(93) 



Cu , Cy — Cu , Xy 



0. 



As we shall see in a moment, the loop equation stays invariant. 

This gauge freedom produces some arbitrary terms in the expectation values of products 
of X operators, which, however, should all cancel out in the momentum superloops and other 
observables, by virtue of periodicity of the physical supermomentum. 

Let us now compute the SUSY area derivative in momentum loop space. The result is a 
SUSY extension of the Mandelstam formula for the gauge potential i X^ 



Ncg'o 
Jo 



2 dC[P]^ 



f exp U J^dPf,iU)xA [X^ , Xy] fexpL dP^{V)X^ 



(94) 



In components it follows from the direct variation. In a superfield language it is even simpler. 
Applying the Qp^(^g^ operator we get sum of terms with 



X, 



Si- 



d6{Si -S) = X^ {5{Si+, -S)- 5{Si^, - S)) . 



(95) 



The first term eliminates the integration over 5"/+!, which yields DPy{S)X^Xy, the second one 
yields —DPy{S)XyX^. The rest of terms add up to the product of two ordered exponentials. 

Now we are in a position to study the superloop equations (|8^). The ^ = equation 
becomes 

r-l 



dp^is) ( 



rexp(i dP^{U)X^ 



+ 



(96) 



Texp U / dP^{V)X^ 



Xu , Xy 



Xx,Xy rexp(z dPf,{W)X^ 







20 



Texp I 



j' dP^{V)X^ 
Jo , 



Texp U / dP^{W)X^ 



Using the factorization of the ordered exponent we can write it as 



dP^(S)rexp ( I I dP^{V)x)j H{X)^^fe^p ( i 



fdP^{W)X, 

J S 



where 



H{X)^^ = + 



{|0) (0| , [Xx,X,]}, 



0, (97) 

(98) 



is some operator, related to X. 

But this is the symmetry property! The transformation 

dX^ = dexHiX)^x, 



(99) 



will result in the same variation of our Ansatz. Therefore, we conclude that the momentum 
loops C [P] are invariant with respect to these transformations with arbitrary infinitesimal 
vector (^e^. Wc could consider the family of operators X^it) parameterized by arbitrary path 
e{t) and obeying the equations 

dtX^, = dtexH{X)^x. (100) 



The correlation functions must be independent of e{t). The derivative 



Nc9o 



dC [P]l 



dt 



£ dP^iS) / fexp L dP^{V)xA H{X)^xdtexfe^p (^t dP^{W)X^ 







(101) 



0, 



for arbitrary non-periodic function P^(S'). 

In the same way, the ^ = 1 equations are equivalent to the symmetry transformations with 



H'{X)^x = S^x + {|0) (0| , [Xx , X,] } [X^ , X, 



The difference of these transformations is the commutator 

H{X)^x - H'{X)^x =[[X^, X,] , {|0) (0| , [Xx , X,] }] • 



(102) 



(103) 



What is the meaning of these observations? The loop equations were derived as Schwinger- 
Dyson equations for the large Nc Yang-Mills theory. Apparently, this theory is equivalent to 
the Quantum Mechanics of the operators X, with some hidden nonlinear symmetry. This is 
not a symmetry of the loop equations, this symmetry is all these equations are about. 
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Any representation of this symmetry solves the loop equations. We could introduce the 
generators T^, T'^ such that 



fA|0)=f'jO)=0, 
(0|fA = (0|f'A = 0. 

The Fock space is defined as the set of words 

\txi...Hm) =al^...al^\Q) , 
(f„...z^i| = {0\at,„...at,^, 

n 

(f„ . . . Z^ll l/il . . . IJ,m) = Smn '^Wr 

i 

The position operator according to Voiculesku can be chosen to be 



(104) 



(105) 



Xfj, — + ^ Qij,,ij,i...ij,i,ali^ ■ ■ -ojifc) 
k=l 



(106) 



where the first term reduces the word by one letter, and the second term expands it by 
arbitrary number of letters. That is 



(f„ . . . ui\ Ifii^i . . . = 1, 

(f„ . . . Vl\X^ It'fc • • • t'n) = Qii,ui...ut,_i, 



(107) 



the rest of matrix elements vanishing. As for the generators, those have arbitrary matrix 
elements, between any finite words. 

Let us write down explicit set of relations between matrix elements. First, consider the 
VEV of these equations, we get 







Xn , Xiy 



Xu , 







(108) 



We used the fact that ( 



X^ , Xu y = 0, which follows from space symmetry in more 
than two dimensions (in two dimensions this could be proportional to e^,^). 

The next equation follows when the matrix elements between vacuum and arbitrary word 
are taken. We get 



(0| X^Tx l/Ui . . 



0. 



All the creation terms in X^ drop here, so we simply find 



(109) 



(110) 
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In other words, not only vacuum state is annihilated by the generator Tx, any (left) letter is 
annihilated as well. For the second generator T'^ we get the similar equation 



t'xX,\0)=0. 



(Ill) 



However, the operator generates the tower of word states from the right vacuum, so that 
this equation is not so simple. 

Finally, the non-vacuum matrix elements of the symmetry relations read 



W 

w 



, Fa 
X^ , F'a 



X^ , Xi, 



0) {[x,,^]\\w') 



(112) 



W) = 6^x6w,W' + {W Xx,X, O) {[f,,u]\\W') 



with obvious notation 



(113) 



Note that these equations are only quadratic in unknown matrix elements of X. All the higher 
order terms canceled. Note also that there is an infinite set of bilinear equations 



W 



Xf_i , Fa 



W 



W 



X^l , T'a 



l-iX"W,W', 



(114) 



for all the words W with n ^ 2 letters. The truly nonlinear equations are those with n = 2 



W 
W 



X^,Tx ap) =6^x{W\\ap) + (W 



X^ , Xp 



X^,T'x ap)=6^x{W\\a(3) + {W Xx,Xp 



5xp 0), (115) 



Xx , Xo 



V/3 



. 



At present we do not know how to solve these equations in closed form. However, the 
recurrent equations for the correlation functions of products of X can be derived and solved 
one after another. 

11 Operator Expansion 

These recurrent equations follow directly from the symmetry property 



^(O X^,...H{X)^^x...X, 



1=1 







0. 



(116) 



At n = 1 we have the previous equation ( [1081 ). In the higher orders we get similar recurrent 
equations. The highest rank tensor here come from the term with / = n. Moving all the rest 
to the right-hand side and skipping the terms which vanish identically, and using the n = 1 
equation we get a recurrent equation 



(0||0)(0 



x^^ . . . ^^„„i 



XiJ.„ , x^ 



Xx , Xy 



(117) 
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1=2 

n-2 

1=3 

n-2 

E 

«=3 



x^^ . . . 




Xp 


X^i . . . 


X^ll 


Xp 







-'^A 5 Xp 







...X, 



...X, 



/in 



The second transformation leads to the double commutator equation 



i:(o 

1=1 



X, 



X^i , Xp 



{|0)(0| , [1a, 1.]} 



..X 



Mil 



) = 0. 



(118) 



There are two highest rank tensor terms, which cancel in virtue of cyclic symmetry. Therefore, 
the second equation reads 

n-2 



Z=3 



X 



Ml • ■ 



X^l , Xy 



,{|0) (0| 



X\ , 



) = 



(119) 



The first terms of this expansion are easy to find. The most general solution for the four 
point function reads 

1 







X^l ^At3 X^^ 







2{d 



symmetric. 



The symmetric terms, proportional to 



(120) 



(121) 



are left ambiguous. These are the gauge terms, which cancel in observables by virtue of 
translational symmetry (p^). The correlation functions of Abelian operators ( C^^ ...C^„) 
produce symmetric terms. 

In four dimensions there is an extra term 

^c5'0^Atl/i2M3/i4 (122) 

in the four point function, coming from chiral anomaly. In two dimension the similar term 



Ncgoe 



A'lM2 



(123) 



would contribute to the two point function. 

The anomalous terms have perturbative normalizations, not affected by effects of strong 
interactions. The normalization can be recovered from (p3[), keeping in mind the extra factor 
of in definition of the vacuum average ( pT]) . 

The two-point correlation in four dimensions remains undetermined. This introduces ar- 
bitrary parameter of dimension of square of length 







X^Xp 







^(0||0) V- 



(124) 
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In the n = 3 order we get the equation for the 6-th rank tensor 











(125) 



The computations of these tensors was done by means of the Mathematica algorithms of 
the previous work (see files attached to the hep-th source of [^). Again, there are some free 
coefficients, which (partly) drop in observables. 

The non-commutative probability theory goes as usual with this definition of X correlators. 
The operator expansion goes as above, with coefficients Q^„^...^^ given by planar connected 
moments ( X,,, ■ ■ ■ X,,„ 



. . . 

I pi. conn 

12 Glueball Equation 

As it was pointed out in |||] the equation for the glueball excitation wave- function ^[C] corre- 
sponds to linearization of the loop equation, which can be achieved by variations of the position 
operator. Let us elaborate this idea in the context of the supermomentum loop equation. 

First, let us identify the total momentum dependence of the glueball loop equation. The 
coordinate superloop equation had the form of 



Lc (W [.]) = Rc (W [.],>V[.]), 

where L is a linear operator, and R is a bilinear one. 
The initial form of the L operator was 



(126) 



Lc (W [.]) 



-0 5X^{0)6X^iS)- 



(127) 



which involved integration over infinitesimal interval (—0, +0) around the point where the 
first variation took place. This operator is local, it picks only the delta terms, related to the 
Y-M equations of motion. 

Now, in virtue of translation invariance 



dS 



5W [C] 
6X^{S) 



0, 



which allowed us to replace the local operator by a nonlocal one 



Lc {W[.]) 



dS 



5W [X]l 



SX^{0)SX^iS)' 



(128) 



(129) 



with integration covering the rest of the loop. This was the operator we used above, in 
transforming superloop equation to the momentum space. 

When the glueball wave function is studied, the translation invariance is replaced by fixed 
momentum condition 



dS 



sxiXs) 



(130) 
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This will result in extra term in the momentum loop operator 



(131) 



Repeating all the previous steps for the linearized momentum loop equation, we find 



rl ( d(c[P]l^[P\\ + ^[P]^C[P]\ 



(132) 



Let us prove that the following Ansatz solves this equation 



B 



B 



dP^^iS) 



TexpU / X,dP,{V)]Y^Texp{t / XxdPx{W) 



provided the operator Yx satisfies the commutation relations 



kxY. 



(133) 



(134) 



The idea is to take a variation 6C [P] = ^[P] of initial operator representation ( pl| ) with 
6X^ = NcQ^Y^, to reproduce the right side of the glueball loop equation ( |132D . The left side 
comes about from the commutator (134). 

The formal proof follows from the identity 







dPf,{S) 



f exp K / X^dP^iV) i;f exp h / XxdPx(.W) , f 



) (135) 



and the same for f^. Commuting the generators with ordered exponentials as follows 



B 



Texp U / X^dP^iV) , r> 



(136) 



dP^iU)Texp I / Xo,dPa{W) H^xiX)Texp i / XpdPp{T) 



B 



it is 



and using the definitions (pq) , ( pl)^ ) , ( |134| ) and algebra of differential forms (pq) , (p 

straightforward to verify the glueball loop equation. 

The linear equations (|134| ) represent the eigenvalue equations for the glueball spectrum in 
this theory. The states 1^ |0) represent the momentum eigenstates, 



f'J'^\Q) = KY^\{}). 



(137) 



We see, that operators fa, f ^ play the role of spatial translation generators of the endpoint 
of "QCD string". They have to be found from above nonlinear operator algebra. In general, 
they do not commute, so these are not the Poincare translation generators. 
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Clearly, the commutators must annihilate these states 



t'^,t'p Y,\0)=0, 
Y,\0)=0. 

which yields restrictions on the Voiculesku coefficients of l^j operators. 



(138) 



13 Quark Propagating Around the Superloop 

Let us now have a look at the negative parity current correlators, which up to normalization 
and standard Chan-Paton factor tr r^^ . . . Ta„ are given by the following integral 



p!^l...!/„ ^^^^ ^ ^ ^ ^^-j 



(139) 



d 











/o 5K,,(5i) 

where it is implied that after taking all functional derivatives the supermomentum 



K^{s) = J2kj,eiSi,s). 



(140) 



i=l 



So, the value of the supermomentum also depends upon Si, . . . Sn- 

There is a systematic way of handling such superpath the ordered exponentials, suggested 



in [11 1 . In our context this looks as follows. We introduce the Quark superfield 

Qw(.S)=qw{s)+eCw(.s) (141) 

which is a \ket) vector in the Fock space of words W = \fii . . . and respectively 

Qw{S)=qw{s) + eCw{s), (142) 

which is a {bra\ vector. Then we represent the SMLoop as follows 







fexp U / dK^,{U)UX^ 







(143) 



The proof essentially repeats the arguments of [|lT|, with some improvements and corrections. 
For brevity we will skip the indexes W, so that, e.g. the propagator 



VQVQexp ( / Q{U) (dQ{U) - idK^{U)X^,Q{U)) j (O |Q(0) Q(l)| 0) . 



G{S,S') = {Q{S)®Q{S'))q^q, 
represents the operator in Fock space with matrix elements 

( W\ G{S, S')\W') = {^ii... ^in\ G{S, s') I/.; 



(144) 



(145) 
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The bare propagator G^{S, S') of this Quark satisfies an equation 

DG^{S,S') = S{S' - S), (146) 

which we know how to solve 

G^{S,S') = e{S,S'). (147) 
The higher vacuum loops all vanish by virtue of the identities 

@{S,S')e{S',S) = 0, (148) 

@{s, s')e{s', s")e{s", s) = o, 

This leaves only one vacuum loop, corresponding to first order diagram 

/' dK^{S)Q{S, S)tTX^ = 11' dK^iS)tiX^ (149) 
Jo ^ Jo 

which vanishes by virtue of periodicity of K^{S), (plus the trace vanishes due to Lorentz 
symmetry) . 

The complete Greens function is a solution of the equation 

DG{S, S') = 5{S' -S)-i {DKf,{S))Xf,G{S, S'). (150) 
The formal solution of this equation is given by the the ordered exponential 

G{S, S') = fexp (^i dK^{U)Xi}j , (151) 

which is just what we need. 

The variations with respect to supermomentum now act on usual rather than the ordered 
exponential with the result 

d {Q{S)Xf,Q{S)) = dQiS)X^Q{S) + QiS)X^dQ{S), (152) 



which plays the role of the vector current of our superloop theory. 
The multicurrent amplitude reads 

C-'^"(fci,...fc„)= (153) 
£ ( {QiS)X{,}Q{S)) ( |Q(0) ® Q(l) I ) ) . 

One may try to integrate by parts here, replacing the covariant derivatives of Quark currents 
by those acting on the G functions in our the ordered integral. However, there is some more 
Si dependence, which comes from external supermomentum K^^S). This brings down more 
Quark currents nQX^Q, which does not simplify formulas. 
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The Gaussian average over Q, Q reduces to super-graphs with propagator G{S, S') in ex- 
ternal field Xfj_DKfj_{S). Now, once all variations are done, we can set the supermomentum to 
its value, after which this external field reduces to the sum of 5 functions in superspace. 

After this the ordered exponential of the integral of sum of the 5 function terms will reduce 
to the product of the vertex operators exp (^i k^^Q{Si)X^Q{Si)^ . 

The flavor vector current correlation functions are given by the multiple superspace in- 
tegrals of the QtQ average. This average reduces to sum of the connected graphs. All the 
vacuum graphs vanish in our ID supergravity, as discussed above. 

The generic connected graph looks the same as the string theory graphs, with currents 
represented by 

r^{k, S)^T''®d {Q{S)X^Q{S)) exp (i k^Q{S)X^Q{S)) , (154) 
and propagators <d{S,S'). 

14 SUSY and Chiral Symmetry 

In the vector sector (anti-periodic boundary conditions) we have to keep the kinetic term in 
the loop Lagrangian. 

The general formula (|4^) for the current correlations involves an extra operator Z, which 
we replaced by 1 in the topological sector. Using the above quark superfield theory we could 
replace this operator by 

Z = lim r dTT"-^ (1 + a In T) exp ( -Tmg - ^ ) , (155) 

L = \£ dSD [q{S)XMS)) {Q{S)X^Q{S)) , 

which will result in extra quartic interaction. 

What could be an effect of this interaction? Formally, one can develop the perturbation 
expansion and compute all the arising T integrals 

r dTT'^-^-' (1 + a InT) exp (-Tml) ^ ( #i±ll _ inml) . (156) 

JO V / n! V ra' / 

We see, that these integrals contain the increasing powers of the bare mass. 

However, the remaining superspace integrals may diverge, because the local SUSY is broken 
down to the global one. Specifically, the singularities at small intervals s — s' ^ now do not 
match the number of integrations. Extra operator = produces ultraviolet divergences. 
Part of these divergences is for real- one has to renormalize the mass after all. 

The spontaneous chiral symmetry breaking must come about as a singularity O 
in normal channels, but also as extra poles at zero momenta in the pion channels of vector 
current amplitudes we consider. 
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This means that the self energy, induced by this quartic interaction, can remain finite in 
the zero mass Hmit, and move the pole to zero momentum. By power counting it is quite 
plausible. The bubble diagrams will contain the singularities, which presumably will be cut 
at the new scale, induced by the mass. Solving the Dyson equation for this new mass scale 
may lead to the chiral symmetry breaking in the limit mo 0. 

The divergences produced by the physical mass term are unpleasant. On the other hand 
we know that the vector sector in the chiral limit is finite. The pion poles involve the finite 
residues which do not need extra renormalization. Is there a way to avoid these spurious 
divergences? 

We suggest the following soft chiral symmetry breaking, without unnecessary violations of 
local SUSY. Let us add the constant Abelian external field with finite density of topological 
charge. The chiral symmetry will be broken by the Adler-Bell-Jackiw anomaly, so that there 
will be no massless pion. 

The vacuum energy can be written as the sum of the zero mode contribution and the finite 
modes. In the limit of the vanishing quark mass 

f = "2 dX iuoS{X) + p{X)) In (mg + A^j (157) 
= —uqIu mo ~ 2 J dXp{X) In yniQ + A j ^ —uq In mo + const, 

where I'o > is the net number of zero modes. The leading divergent term in the chiral limit is 
the contribution of the zero modes. Unlike the chiral anomaly, the net number of zero modes 
is not a topological invariant. 

Let us study this number in the free quark theory. The corresponding Euler-Heisenberg 



Lagrangian was computed in Appendix B. The formula (170) can be written as 

V drr"-i(l + alnT)e-^™o]^6.coth(26iT), (158) 

i=i 

where =bz6j are the eigenvalues of B^,^. In the chiral limit the leading logarithmic term comes 
from the region b^^ <C T <C ttiq"^ where one could replace coth(26jT) — > sign{bi), after which 
the integral yields 



Jo 



dTT^-i (1 + a In r)e-^""o ^ _(^ + In mg)U VdetzB . (159) 

This is an absolute value of the Pfaffian. 

The same result (up to normalization) could be obtained much simpler, by neglecting the 
kinetic energy in the path integral. The path integral becomes super symmetric, 

(160) 



V I VX5^{X{0))exp i^-B^, j\x^{S)X,{S) 
V j PX5''(X(0))exp VY.h j dX2i^i{S)X2i{S) ocUjJ|6i| = v\\/deiiB 
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We used the rotational invariance to reduce the antisymmetric quadratic form to the Jordan 
form. Then we rescaled half of components of the superfield X2i{S) —>■ b~^X2i{S). The 
SUSY measure stays invariant, because STrl = 0, therefore the only contribution comes from 
rescaling of the S function. The proper time integral gives then the logarithmic contribution. 

So, we get the absolute value of the topological charge density. The vacuum average of this 
number in QCD does not vanish, unlike the average topological charge. In the background of 
external constant flavor field with nontrivial topology wc expect this number to be proportional 
to the volume of space, which makes it a leading term in the chiral limit (it is proportional to 
Inrng, unlike the contribution from the finite modes). 

We could argue that the quark kinetic energy drops in QCD by the same reason as in the 
Euler-Heisenberg Lagrangian, and we arrive at the locally SUSY theory in the vector sector 

r;;i---"(fei,...fe„)= (161) 

-In mo [{ (0\Q{0)^Q{1)\0) f\d{Q{Si)X,,Q{Si))) . 
with the usual supermomentum 

K^{U) = j2k^^e{Si,U). (162) 

i=l 

The corresponding term in the Quark action also represents the quartic interaction 



B^, J dSD (^Q{S)X,Q{S)) {Q{S)Xi,QiS)) , (163) 



which preserves the local SUSY. The only symmetry which gets broken is the Lorentz sym- 
metry, as we now have external tensor field. 

One could study the perturbation expansion in this quartic interaction and try to identify 
and sum up the terms responsible for the spontaneous chiral symmetry breaking in the limit 
0. The good thing about this perturbation expansion is that all the graphs will still 
be locally SUSY, and therefore calculable. 

This interesting issue calls for a further study. 

15 Summary and Conclusions 



Let us summarize the status of the superloop theory of the large Nc QCD. Here are the 
basic features of this theory. 

• The chiral momentum loop amplitudes, with projector ^ (1 ± 75) inside the Dirac trace, 
satisfy the closed equation, generalizing the ordinary momentum loop equation, for the 
case of the supermomentum loops. 
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The supermomentum loop equation can be interpreted as a nonlinear symmetry relations 
( [L04D of the position operator of the QCD string. This symmetry lead to recurrent 
relations for the coefficients of the Taylor expansion of momentum superloop. 



• The (non-commutative) generators f^,f^ of these symmetry transformations can be in- 
terpreted as string endpoint translation operators. In particular, the glueball excitations 
with fixed momentum are related to the eigenstates of these operators with common 
vector eigenvalue k^. 



ground of strong constant electro-magnetic field are related to (functional derivatives 



• One could introduce the Quark superfield at the loop Q{S) in such a way that the fiavor 



• In the chiral limit the Quark superfield is Gaussian, with the propagator, depending on 

and momenta of all the currents. The exact operator solution for this propagator 
looks very similar to the good old string diagrams. 

• The low-energy expansion of fiavor current correlators in the chiral limit in the back- 
ground of strong constant electro-magnetic field is also calculable in terms of these 
coefficients. There are no singular integrals involved. All the superloop integrals reduce 
to rational numbers. These tedious calculations can be simplified by means of generating 
function derived in Appendix E. 

• In order to account for the chiral symmetry breaking effects in anomalous channels of 
these current correlators, one should add the quark mass niQ. This results in extra quartic 
interaction in our ID superfield theory. The interaction constant is proportional to the 
m-Q. This interaction breaks the local SUSY down to the global one, and introduces the 
extra divergences. Presumably, summing the leading divergences (bubble graphs?) will 
account for the chiral symmetry breaking. 
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• The correlators of the vector quark currents ( Yl J^ii^^) ) the chiral limit in the back- 



of) the SMLoop with supermomentum K^{S) = J2i ^^©(-S'i, S). 



vector current is represented at the loop as ( |154| ). 
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A One More Representation of a Log 



Let us consider the identity 

A-"r{a) = / dTT'^-^e-^^, (164) 
Jo 

valid for any positive definite operator A. Multiplying it by a ant taking derivative ^ we find 

;>oo 

^-°r(l + a)(V^(l + a)-ln^) = / dTT^-^ (1 + a In T) e"^^. (165) 

Jo 

Tending a — > +0, 

7 + lnA = -lim / dTT"'^ (1 + a In T) e"^^. (166) 

a^+OJo 

The formula in the text corresponds to j4 = mo + (z7/iV^)^ which is positive definite regardless 
of the possible zero modes of the Dirac operator. 

B Euler-Heisenberg Lagrangian 

Let us compute the vacuum energy in constant Abelian field B^{X) = | B^yX^. 

The quadratic forms in exponential of (p^) are the same for Bose and Fermi parts 



C ds {x'^{s)Q^, (^^) x,{s) + V/.(s)Qa<. (^) Ms)) , (167) 



with the only difference coming from the boundary conditions]^ 
Let us expand the fields in Fourier expansion on the circle 



X 



(s) = ^ x)] exp (2m vrs) , (168) 

— oo 
oo 

V'/.(s) = 5I^Mexp((2n + l)z7rs), 

— oo 

then the functional integral yields the product of ratio of the determinants for each Fourier 
mode 

oo VdetQ(2m7r) 

The first factor corrects the contribution from the zero mode x^. There is, in fact, no 
Gaussian integration for this mode in our action, so that this integration produces the volume 
V of space. The factor ^det Q{0) removes the zero mode factor from the denominator at 
n = in the product. 



■^Let us note in passing that the periodic boundary conditions V'(l) = V'(O) correspond to the chiral anomaly 
instead of the vacuum energy. We discuss this issue later in some detail. 
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Now, the standard computation of the infinite product yields 

V dTT'^-^ (1 + a InT) e-^^'o-^det {B cot 2BT). (170) 

which is our representation of the Euler-Heisenberg Lagrangian. The determinant of the 
matrix function is imphed here. One could expand it in power series in B(we denote the 
Bernoulli numbers as B2n) 

Boot 2BT = ^ (^1 - E |4|(4^^)'"J ' (171) 



^det {B cot 2BT) = T'^ ^ C2n{B)T^'' 

n=0 



with some invariant polynomials 

C2n{B) = J2mI[(^^B'T- (172) 

{k} 

The proper time integrals 

dTT"-^ (1 + alnT) g-^'^oT^"-!'^ = ^ (^amJ^-^^-^^r (^2n + a - ^dj^ . (173) 

The pole at d = 4 in the B^ term is the usual UV divergence, removed by our a renormalization 
as follows 

— (amo (a)) ^ -7 - 2 In mo- (174) 

C Supersymmetrization Theorem 

Let us consider the function F {Si, . . . Sn) of n variables Si = (sj, 9i) which linearly transforms 
under infinitesimal local SUSY transformations 

ds = eep{s), (175) 
69 = eP{s). 

The finite transformation have the form 

s{\) = s + \9P{s), (176) 

9(X) = 9 + A/3(s) + ^X^P{s)p'is). 

One may verify that these polynomials satisfy differential equation 

s'{X) = 9{X)(3isiX)), (177) 
9'{X)=P{s{X)). 
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The local SUSY implies that the equation F {Si{\), . . . Sn{X)) = is, in fact, independent 
of A for arbitrary Grassmann function /3(s). Using this independence, we may choose /?(s) so 
that it "passes through all the 9 points" 

I3{si) = Oi. (178) 
An explicit example of such function is a Lagrange interpolating polynomial 

Now, taking A = — 1 we observe that 

ei{-l) = 0,Si{-l)=Si. (180) 

So, the function, which vanishes identically at the real axis, will vanish identically in whole 
superspace. Therefore, one can extend in superspace any invariant identity which holds on a 
real axis. 

Being valid at any finite n, in a limit it must be also valid for ID functional F[r] which 
depend on some curve F : S = {s{t),9{t)) in superspace. The same transformation with 
A = — 1 moves this curve to the real ax;is provided 

Pis{T)) = e{T). (181) 

Any solution of for inverse function r(s) will give an example of such transformation 

m = dir{s)). (182) 

Such an obvious theorem must, of course, be well-known to mathematicians, but for a lazy 
physicist it was easier to "prove" it than to dig into the mathematical literature. 



D Factorization of Fourier Integral 

The factorization of the Fourier exponentials 

exp (^i dP^{U)X^{U)^ = exp (^i j\p^{y)X^{y)^ exp j\p^{W)X^{W)^ , 

(183) 

follows from the additivity of the dX^ form. When we compute the gradient precisely in a 
breaking point, in a complete integral it is done by means of the integration by parts, with 
the usual result 

exp (^-i 1^ dP^{U)X^{U)^ exp dP^{U)X^{U)^ = -idX^{S). (184) 
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DX^iS) = -DX,{S-) + -DX,iS+). (186) 



In a split form one can obtain the same result in one plays by the rules. This computation 
goes as follows 

-^ exp (^-^ dP^{V)X,{V)^ exp (^i j' dP^{V)X^{V)^ = (185) 

^ d5{V - S)X^iV) = 6iO)X^{S) - S{A - S)X^{A) - dX^{V)S{V - S) = 

A J A 

- dX^{V)6{V -S) = - J dVDX^{V)6{V - S)e{A, V)e{V, B) = 
-DX^{S)Q{A,S)Q{S,S) = -]^DX^{S)Q{A,S). 

The second exponential gives another half of the derivative DX^{S). 

The same result could be obtained much simpler, using the Fourier ^ rule 

^DX^iS-) + \. 

This ^ rule is a general rule to interpret the discontinuities. 

The terms with S— {S+) is the boundary value of the ordinary derivative, acting on the 

left (right) exponential, where they can be replaced by functional derivatives ^^^-^ for the 

inside points. The left derivative |q_>^ is equivalent to j^j^, so it acts on C [P]\, whereas 

the right one acts on C [P]q . In the same way, the left derivative ^{^s-) ^^^^ ^ [-^lo' 

the right one acts on C [P^g- As a result we get all four possible terms, each with the weight 
1 

4- 

E Generating Function for Super- Graphs 

There is a simple solvable model which serve as generating function for these super-graphs. 
Consider the sum of two Abelian constant gauge field strengths G^i, + B^i, . The superpath 
integral is computed in Appendix B, and it reduces to a Pfaffian 

J VX6''{X{0)) exp Q {G^, + B^,) dX,{S)X^{S)^ oc ^det(i? + G). (187) 

On the other hand, we could replace the G simplectic form by an oscillator commutation 
relation 

[R^,R,] = iG^,, (188) 
and compute the vacuum average instead 

' " ^ ^ 'o\ (189) 

/ R 



J VXd'^iXiO)) ^0 f exp (^J^ dXf,{U)R 
exp (^-B^, j\x,{S)X^{S)j oc ^dei{B + G). 



This formula (see Q) follows from the definition of the the ordered product and the Baker- 
Hausdorff formula. Another proof is to switch to Schwinger gauge in the ordered exponential. 
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treating i?^ as a non-Abelian gauge field R^{S). This field is constant, but it has a finite field 
strength ( 188 ). This field strength being Abelian, the Taylor expansion in a Schwinger gauge 
terminates, so that one may replace 



R,iS)^ -G^^X^iS), (190) 



which recovers the initial Abelian formula. 

Comparing the terms of expansion in powers of G we get relations 



{0\R^,...R,,JO) f {<f^X^^}{{U})) = (191) 

Jq \ ^ ' I Gauss 



E^^tr(5-G)^ 



exp 



with the Gauss averages 



{DX^iS)DX,iS') ) = tB^^D6{S - S'). (192) 



On the right-hand side the term O (G") is left in the Taylor expansion. The vacuum average 
can be computed from commutation relations, which gives the set of sum rules for the locally 
SUSY integrals we need. 

The point is this set of relations is valid for arbitrary d, while the SUSY integrals are 
universal numbers, independent of d. So, these relations do not terminate at n = 4 in four 
dimensions, where we need them. 

Let us be more specific. Expanding the exponential on the right side, we get 

E E • • • E *.E'. n ' ■ (193) 

1=0 ''■ kx=l h=l i=l 



Each factor of corresponds to the Gauss contraction (192), and each factor of G corre- 
sponds to the commutator (|18^ ). The indexes between the and G tensors are contracted 
according to the specific Gauss diagram. The coefficient in front of this product of tensors 
B-^ (g) . . . B-^ (g) G (2) ... G on the left is given by the SUSY integral 

± j\'^''UY[D5iUi-Uj). (194) 



. f (-1)"'' ni _ 



The same coefficient on the right is given by above sum of — — nj=i 2F' independently of 
d. In particular, the first relation reads 

-ivB-^G = -iB-^ {Q\RMQ) / d^UD5{Ui - U2) = (195) 
2 Jo 

-^z B-^ {0\[R^, R,]\0) d^UD6{Ui - U2) = 

]:B-^G^u C d^UD5{Ui - U2) = -UtB'^G d^UD5{Ui - U2), 
^ Jo ^ Jo 
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from which we get 

/ (fUDS{Ui - U2) = -1. (196) 

Jo 

These relations in higher orders leave some undetermined SUSY integrals. 

The typical integral involves the products of the O functions and their derivatives, i.e. 5 
functions. The number of S functions is twice less than number 2n of integration variables, 
moreover, there are 6 functions for every variable, otherwise the local SUSY would be violated. 

The simplest graphs are those, when each 6 function is differentiated only once. In this 
case all the points will collapse to one of the endpoints of the interval 0, 1. Some @ function 
will collapse as well, leading to a bubble graph at one or another endpoint, 



but some other will connect the endpoints, contributing just 0(0, 1) = 1. Each bubble 
will give 6(0,0) = 0(1, 1) = ^, so that the integral would be ±2~^ where k is the number of 
bubbles. 

The generic term which arises in momentum expansion involves derivatives of the 5 func- 
tions 

± £ (PS £ <rK^^} {{U})Y[ DS{Ub - Sc) (197) 

where {1 . . . m} is partitioned into {a}, {b} and {c} = {1, . . . n}. There are precisely n indexes 
in the {b} list, and m — n in the {a} list. 

Substituting here the above definition of K^{U) we get superposition of terms with prod- 
ucts of Cronecker deltas times 

£ (PS £ (TU n 5{Ua - Si) n D5{Ub - Sc) = (198) 

/OO poo poo poo 

dSi... dSn dUi... dUm 

-00 J— 00 J— 00 J— 00 

9(0, Si)... eiSn, 1)6(0, Ui)... e{Um, i) H ^(^^ - s,) n mub - s^) 



These superspace integrals are universal numbers, which depend on the partition of indexes. 
This follows from the local SUSY, which is preserved here. 

The derivatives D in the last line can be integrated by parts, cither over Sc or over Uh in 
such a way, that they act only on the 6 functions, rather than 6 functions. We get then 

/oo /"OO POO roo 

dSi... dSn dUi... dUm n Q{Su 5hi) n 0(^ik' Uk+i) (199) 
-00 J —00 J —00 J —00 1 
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n mj - Uj+i) n ^{Sr - Sr+i) n ^(Ua - Si) n m - sc) 

The number of 5 functions equals m + n, therefore all the integrals can be done. In case 
there are cycles, there will be some §{0) factors. All such terms drop, as 

(5(0) = 0, (200) 

S{S - S')S{S' -S) = 0, 

5{S - S')6{S' - S")6{S" - 5) = 0, . . . 

in superspace. So do the terms with the cycle of @ functions, as we already discussed above. 

These identities can be proven by explicit calculation. 

So, only terms without cycles will be left. In those terms we shall have the products like 
d{Ui — C/j_|_i)6(?7j, f/j_|_i). These products can be replaced by ^5{Ui — C/j+i), since &{U, U) = ^. 
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